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ABSTRACT 

A Dirac particle is represented by a unitarily evolving state vector in a Hilbert space which 
factors as H sp i n ® imposition- Motivated by the similarity to simple models of decoherence 
consisting of a two state system coupled to an environment, we investigate the occurence of 
decoherence in the Dirac equation upon tracing over position. We conclude that the physics 
of this mathematically exact model for decoherence is closely related to Zitterbewegung. 
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Decoherence — the process by which a quantum system initially in a pure state evolves 
into a mixed state as it interacts with its environment — has been studied both in toy models 
designed to elucidate aspects of quantum measurement [1-4] and in more complex ones 
modelling realistic physical situations [3,5-7]. In the purely quantum mechanical models 
the Hilbert space factors as a tensor product H = H sys <g> H env of Hilbert spaces describing 
the degrees of freedom of the system and environment, respectively. From the complete 
state \*f>) G H one forms the density operator p := ® which is traced over H env 
to give the reduced density operator p := Tr env p. Physically, the trace corresponds to 
the absence of measurements on the environment; p is the average over all possible states 
thereof. Generically, p is not of the form |\I/ sys ) <g> (* S ys| for any \^ sys ) G iisys; it is rather 
a linear combination YliPi^i) ® (V^l) where G H sys and the coefficients pi are positive 
and sum to 1. That is, p describes a mixed rather than a pure state, the degree of mixing 
being measured by the entropy S := — Trplogp. 

The simplest decoherence models considered typically have a two dimensional H sys 
coupled to an environment whose Hilbert space is also two [1-3] or at most finite [2,4] 
dimensional, although more realistic (infinite dimensional) environments have been used 
to investigate the absence of chiral superpositions in molecules [3], for example. Now, 
the Dirac equation describes the unitary evolution of a spin--| particle — the Hilbert space 
factors as H = H sp i n ® imposition — and H sp i n is finite dimensional. In the spinor chain path 
integral formulation of the Dirac equation [8,9], to each direction in space is associated a 
spin projection, and over a time interval 5 a projection operator is understood to propagate 
the particle a distance 5 in the associated direction. A single path contributing to a 
transition amplitude from x to x' consists of a sequence of spin projections and chiralities 
such that the sum of the translations associated with the spin projections, each multiplied 
by the corresponding chirality, is proportional to the difference in position x' — x, where 
the proportionality constant depends on the spatial dimensionality. The amplitude of such 
a path is proportional to (idm) r (in units with h = 1 = c), where m is the mass of the 
particle and r is the number of chirality reversals. 

From the point of view of decoherence models, this path integral describes exactly how 
the spin degrees of freedom are 'recorded' by the position of the particle. That is, imposition 
can be thought of as the Hilbert space of an environment measuring a quantum system with 
Hilbert space H sp [ n . Thus the Dirac equation can be analyzed for decoherence via reduced 
density matrices. In 1 + 1 dimensions* the particle wave function is a two component 
spinor field ip(x,t) = Oip-i(x,t),^+i(x,t)) in the chiral basis, so the density function is 

p(x, a, x', a', t) = ip a (x, t)ip a >(x', t), for a, a' G {±1}. Tracing over the environment is here 
an integration over position, so the reduced density matrix has four components 



In this note we examine (1) for decoherence by computing, exactly, its entropy as a function 
of time. 

* Since there is no spin — only chirality — in 1 + 1 dimensions, this will simplify the subsequent discussion 
with no conceptual loss and allow direct comparison with decoherence of other two state systems [1-4] . 
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Using the path integral formulation of the Dirac equation it is straightforward to 
compute the propagator K a2Ctl (x 2 , t 2 , t\) satisfying 

rx 2 +At 

^ a2 (x 2 ,t 2 ) = K ct2Ctl (x 2 ,t 2 ,x 1 ,t 1 )yj ai (x 1 ,t 1 )dx 1 , (2) 

Jx 2 -At 

where At := t 2 - h. The result is [8,9,10]: 

„ , . \ f S(At — aiAx) — (At + a\Ax)mJ\(mr) jlr if a 2 = cki; ro . 

Ka ^ X2 ^ XlM) -\imJ Q {mr)/2 if a 2 ± a u (3) 

where Ax := x 2 — x\, r := a/ (At) 2 — (Ax) 2 and Ji denotes the i th order Bessel function 
of the first kind. Using (2) and (3) in (1) gives an exact formula for the reduced density 
matrix as a function of time in terms of the initial state of the particle. 

Decoherence is usually studied for an 
initial state which is a tensor product pure 
state l^sys) <8> |^env) [1~7], so consider the 
case ip(x,0) = where f(x) is pro- 

portional to e~ x2 / 2 and is normalized so 
that fi/>*(x,0)if>(x,0)dx = 1. That is, 
the particle is initially in an equal super- 
position of negative and positive chiral- 
ities and is localized near x = 0. Fig- 
ure 1 shows the entropy S(t) of the re- 
duced density matrix (1) for the range of 
masses m G {0, 1, 2}. Near t = the rate 
of entropy gain increases as the mass de- 
creases: the top curve is for m = while 
the bottom one is for m = 2. Thus, as 
could have been anticipated from the am- 
plitudes associated to paths in the spinor chain path integral, m _1 controls the strength of 
the system/environment interaction: when m _1 is larger, decoherence occurs more rapidly. 

To understand this decoherence, we begin by considering m = 0. In this case the Bessel 
function terms in (3) vanish so the effect of the propagator is simply to translate the initial 
chirality position amplitudes along the lightcone: ip(x, t) = f(x + t) (J) + f(x — t) (°). For 
the initial position Gaussian we can compute the reduced density matrix (1) exactly to 
get: 

« t) = Ke-' a e ~0' <4) 

The diagonal terms in (4) are constant and equal, while the off-diagonal terms go rapidly 
to 0; thus the entropy goes rapidly to 1 as indicated by Figure 1. 

Next, we consider m = 1. Just as in the massless case, if we were to plot the position 
distributions for each chirality at t = 0, i.e., \ip-i(x, 0) | 2 and \ip + i(x, 0)| 2 , they would 
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Figure 1. Entropy as a function of time for the 
reduced density matrix (1) with the particle in 
an initial equal superposition of chiralities and 
localized in a Gaussian around x = 0. From top 
to bottom, the curves are for increasing mass. 
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Figure 2. The t = 1 chirality position distribu- 
tions for m = 0. The distribution with the peak 
at the left/right is for negative/positive chirality; 
each remains a Gaussian. 



Figure 3. The same chirality position distribu- 
tions shown in Figure 2 but for the case m = 1. 
Dispersion distorts the initial Gaussian and re- 
duces the separation of the distributions. 



be two coincident Gaussians centered at x = 0; there is no initial correlation between 
position and chirality for a tensor product state. Figures 2 and 3 show the chirality 
position distributions at t = 1 for m = and m = 1, respectively. In each case the 
distributions separate: a particle with negative/positive chirality tends to move to the 
left/right. Now chirality is correlated with position, more strongly for m = than for 
m = 1, and it is this correlation which is measured by the entropy upon tracing over 
position. For m = 1 dispersion distorts the initial Gaussian so that there is a weaker 
correlation between chirality and position than for m = 0, and thus lower entropy. 

Since the spin system is coupled to an infinite dimensional position Hilbert space 
environment, not a small finite dimensional Hilbert space with a short Poincare recurrence 
time [2-4], one might have expected even this weaker chirality/position correlation to 
continue to build and the entropy to increase asymptotically to 1 just as it does for m = 0. 
This expectation is contradicted, however, by the m = 2 curve in Figure 1. In this case the 
entropy reaches a (local) maximum and then begins to decrease. In fact, the oscillatory 
behavior beginning here characterizes all of the massive entropy curves. 

To explain how the entropy can decrease, despite the infinite dimensional position 
Hilbert space, Figure 4 shows the m = 1 entropy curve over a longer time interval < t < 2. 
The inset graphs show the same chirality position distributions plotted in Figures 2 and 3, 
for times t e {0.5, 1, 1.5, 2}. The slowing and then reversal of the entropy increase is due 
to the decreasing portion of the distributions in the fastest moving (outermost) peaks and 
the resulting increasing portion of each distribution which overlaps the other.* That is, 
dispersion acts initially to increase (by separating the distributions) and subsequently to 
decrease (by broadening/distorting the distributions) entropy; then the cycle repeats. For 
smaller masses not only does the initial increase occur faster, but the entropy also reaches 
greater values before starting to decrease. 



* Spacetime plots of quantum lattice gas automaton simulations of the Dirac equation showing exactly 
this behavior (although with periodic spatial coordinate) can be found in [10]. 
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Figure 4. Entropy plotted as a function of time for an initial position Gaussian with m = 1. 
The inset graphs show the chirality position distributions at half integer times. 

The results shown in Figures 1-4 have been obtained for quite special initial con- 
ditions: in each case an equal superposition of negative and positive chiralities. While 
this choice demonstrates decoherence in the Dirac equation particularly clearly, any spinor 
tensored with an initial position Gaussian will decohere. Figure 5 shows the entropy curve 
over < t < 1 for initial condition yj(x,0) oc e~ x ^(J with m = 1, while Figure 6 shows 
the chirality position distributions at t = 0.5. The larger distribution is \ip + i(x, 0.5)| 2 ; as 
expected, since the initial chirality is purely positive, the chirality is still more likely to be 
positive than negative at t = 0.5. Nevertheless, there is some probability that the particle 





Figure 5. Entropy as a function of time for 
the chiral initial condition localized by a position 
Gaussian and for m = 1. Initial entropy gain is 
slower than in the corresponding parity invariant 
case shown as the middle curve in Figure 1. 



Figure 6. The distributions |V>±1 (x, 0.5) | 2 for 
the same initial condition as in Figure 5. The 
larger distribution is for positive chirality; the 
smaller, for negative. There is a small chiral- 
ity/position correlation. 
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has negative chirality, so chirality /position correlation has developed. Just as in the parity 
invariant cases, only more slowly, this correlation causes the entropy increase shown in 
Figure 5. 

By superposition in H sys , therefore, any tensor product initial condition of the form 
ip(x,0) oc e~ x / 2 (a,6) decoheres. This result, and its familiarity in the context of other 
models of decoherence [1-7], might lead us to believe that any tensor product state, initially 
localized in position, decoheres — and to conclude that the spin of a Dirac particle evolves, 
absent position measurement, stochastically rather than purely quantum mechanically. 
This, of course, is not true, as we shall see by working in the stationary state basis. 

The stationary states of the Dirac equation are the plane waves ip( k > e \x) = e lkx u(k, e), 
where k G K, e G {±1}, and the fixed spinors u(k, e) G C 2 are normalized and orthogonal 
for e = 1 and —1. While it is not localized in space, of course, suppose we take a plane 
wave as our initial tensor product state. Since this stationary state evolves by 

multiplication by the phase e~ l£Ujt , where the energy ecu satisfies the dispersion relation 
u 2 = k 2 + m 2 , the density operator is constant: 

p(x, x', t) = e lk{x - x,) u(k, e) <g) u\k, e) =: e ^ x ~ x ') p{t), 

as is the reduced density matrix. Thus stationary states do not decohere, a familiar result 
in decoherence models with tensor product eigenstates of the Hamiltonian [1,3,4]. 

Any initial state can be expanded as a superposition of plane waves, i.e., 

i;(x,0) = J2 [ dk^(k,e)e ikx u(k,e). 
In this form the time evolution is transparent: 

ip(x,t) = ^2 f dkyl{k,e)e lkx - ieujt u{k,e), 

so it is easy to write down the time dependent density operator explicitly: 



To obtain the time dependent reduced density matrix we need only set x' = x in (5) and 
integrate over x. The integral of e^ h ~ h ) x gives a delta function of k — k', so we can evaluate 
the integral over k' to get 



p(t) =J2j dk^(k 7 e)e- l£Ujt u(k,e)(g)u^(k 1 e')e ie ' ujt ij(k 1 e'). (6) 

e,e' 
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The only nontrivial time dependence of p(t) derives from terms in (6) with a factor 



is, decoherence depends on the presence of positive and negative energy modes of the 
same momentum. The initial states whose decoherence we exhibited in Figures 1-5 are 
superpositions of such modes; in contrast, a wave packet constructed from only positive 
energy plane waves, would disperse, of course, but not decohere. 

We conclude by remarking that while the Dirac equation provides a mathematically 
exact model of decoherence, this decoherence is another aspect of the physical difficulties 
with the one particle interpretation of the Dirac equation. The particle decoheres and be- 
haves stochastically rather than quantum mechanically exactly to the extent that it is not 
identified as an electron or positron, for example. The same interference between positive 
and negative energy states leads to Zitterbewegung [11]. That the entropy oscillates, rather 
than increasing permanently as it does in other decoherence models with infinite dimen- 
sional environment Hilbert spaces [5-7], is a consequence of the same oscillating factors 
e -2ieujt j n ^ £ na j. cause Zitterbewegung. To revise FitzGerald's rendering of The Rubaiyat 
of Omar Khayyam [12]: 




The Moving Finger writes; and, having writ, 
Zitterbewegt sich. 
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